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Abstract 

Within the QCD sum-rule approach in an external field, we calculate the 
baryon matrix element of isovector-scalar current, H-q = {B\uu — dd\B)/2M-Q, 
for octet baryons, which appears in the response of the correlator of baryon 
interpolating fields to a constant isovector-scalar external field. The sum 
rules are obtained for a general baryon interpolating field with an appropri- 
ate form for the phenomenological ansatz of the spectral density. The key 
phenomenological input is the response of the quark condensates to the ex- 
ternal field. To first order in the quark mass difference 8m = m^ — m u , the 
non-electromagnetic part of the baryon isospin mass splitting is given by the 
product of 5m and Hb- Therefore, QCD sum-rule calculation of Hb leads to 
an estimate of the octet baryon isospin mass splittings. The resulting values 
are comparable to the experimental values; however, the sum-rule predictions 
for Hb are sensitive to the values of the response of the quark condensates to 
the external source, which are not well determined. 
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I. INTRODUCTION 



To understand the observed properties of hadrons from the underlying theory of the 
strong interaction, quantum chromodynamics (QCD), is a challenging task since the QCD 
remains intractable at low energies. Among the attempts made in dealing with the strong 
interactions at low energy scales is the QCD sum- rule approach fll], which has proven to be a 
useful tool of extracting qualitative and quantitative information about hadronic properties 

HI- 

One of the extension of the sum-rule methods made by Ioffe and Smilga for external 
field problems enables one to calculate the baryon matrix elements of various bilinear quark 
operators. These include the matrix element of electromagnetic current to determine the 
magnetic moments , the matrix element of the axial vector current to find the renormal- 
ization of baryon axial coupling constant ||[7|], the matrix element of the quark part of the 
energy momentum tensor, which gives the momentum fraction carried by the up and down 
quarks in deep inelastic scattering [§,0, and the matrix element of isoscalar-scalar current 
for evaluating the nucleon sigma term |10| . 

In this paper, we evaluate the baryon matrix elements of isovector-scalar current, 



Hb = (B\uu — dd\B) /2Mb, within the external field QCD sum-rule approach. In Ref. |TT 



the proton matrix element (p\uu — dd\p)/2M p has been calculated in the external field ap- 
proach. However, one piece of the phenomenological representation has been omitted in the 
calculation. This has been pointed out recently by Ioffe |12[. In the present paper, we shall 
derive the appropriate phenomenological representation, which includes the piece neglected 
in Ref. We use this complete phenomenological representation and a general baryon 



interpolating field to calculate the matrix element H-q for octet baryons. 

External field sum rules for baryons are based on the study of the correlation function 
of the baryon interpolating field in the presence of an external field. The appearance of 
the external field leads to specific new features in QCD sum rules which distinguish them 
from those in the absence of the external field. At the hadron level, the spectral parameters 
usually used in the parametrization of the spectral density, baryon masses, pole residues, and 
continuum thresholds, all respond to the external field. Consequently the phenomenological 
representation for the response of the correlation function contains a double pole at the 
baryon mass whose residue contains the matrix element of interest. This corresponds to 
the response of the pole position. The response of the pole residues gives rise to single 
pole terms, which contain the information about the transition between the ground state 
baryon and excited states. The single pole contributions are not exponentially damped after 
Borel transformation relative to the double pole term and should be retained in a consistent 
analysis of the sum rules. In addition, there are terms corresponding to the response of the 
continuum thresholds, which should also be included in the calculation. At the quark level, 
the external field contributes in two different ways-by directly coupling to the quark fields 
in the baryon current and by polarizing the QCD vacuum. By equating these two different 
representations for the response of the baryon correlator, one obtains the external field 
sum rules, which relate the baryon matrix elements of various current to QCD Lagrangian 
parameters, vacuum condensates, and the response of condensates to the external source. 

The observed baryon isospin mass splitting has its origin in the electromagnetic interac- 
tions between quarks and in the different masses of the up and down current quarks. The 
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contributions of the latter to the first order in the quark mass difference 5m = m^ — m u is 
given by the product of 5m and the baryon matrix element of the isovector-scalar current. 
Therefore, the QCD sum-rule calculation of the baryon matrix elements of isovector-scalar 
current for octet baryons naturally leads to an estimate of the octet baryon isospin mass 
splittings. (The S _ — S° and S + — S° splittings will not be considered here as there is 
mixing of the E° with the A via isospin- violating interactions.) 

The rest of this paper is organized as follows. In Sec. [TT], we establish the baryon QCD 
sum rules in an external isovector-scalar field. In Sec. |HI| , we then analyze the sum rules 
and present the results. In Sec. IV], we estimate the isospin mass splittings of the octet 
baryons using the baryon matrix elements of isovector-scalar current calculated from QCD 
sum rules. Further discussion of our results are given in Sec. [V] 

II. BARYON QCD SUM RULES IN AN EXTERNAL ISOVECTOR-SCALAR 

FIELD 

In this section, we establish the baryon QCD sum rules in the presence of an external 
isovector-scalar field. In previous works ]3]-|T2"||, the phenomenological representation for the 
correlator is usually obtained by analyzing a double dispersion relation. Here we present 
an alternative approach to derive the phenomenological representation. The operator prod- 
uct expansion (OPE) results can be easily obtained following the procedures outlined in 
Ref . []l l]] . We work to leading order in perturbation theory and to first order in the strange 
quark mass. Contributions proportional to the up and down current quark masses and the 
gluon condensate are neglected as they give numerically small contributions. We include 
condensates up to dimension eight. 

Consider the correlator of the baryon interpolating field in the presence of a constant 
external isovector-scalar field Sy 

n(SV,g) = i J d 4 xe^(0\T[ VB (x)rJ B (0)}} Sv , (2.1) 

where i]b is the interpolating field for the baryon under consideration. We consider baryon 
interpolating fields (currents) that contain no derivatives and couple to spin-| states only. 
There are two linearly independent fields with these features, corresponding to a scalar or 
pseudoscalar diquark coupled to a quark. In this paper, we take a linear combination of 
these two fields, 

r] p (x) = 2e abc \[u^{x)Cd b {x)\^ h u c {x) + t[u^(x)C^ 5 d b (x)]u c (x)^ , (2.2) 
r]Y,+ (x) = 2e abc [[u T a (x)Cs b (x)]-i 5 u c (x) + t[u^(x)C-f 5 s b (x)]u c (x)j , (2.3) 
T] s o(x) = 2e abc {[Sa(x)Cu b (x)]>y 5 s c (x) + t[s T a (x)C^u b (x)}s c (x)} , (2.4) 

where u(x), d(x) and s(x) stand for the up, down and strange quark fields, a, b and c are 
the color indices, C = —C T is the charge conjugation matrix, and t is an arbitrary real 
parameter. The interpolating fields for neutron, S~ and S~ can be obtained by changing 
u(d) into d(u). The interpolating fields with t = —1, advocated by Ioffe |T3"|,|ri|] , have 
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been used exclusively in previous papers on external field sum rules 0-0- 111 principle, 
the sum-rule predictions are independent of the choice of t; in practice, however, the OPE 
is truncated and the phenomenological description is represented roughly. The goals in 
choosing the interpolating field for QCD sum-rule applications are to maximize the coupling 
of the interpolating field to the state of interest relative to other (continuum) states, while 
minimizing the contributions of higher-order terms in the OPE. These goals cannot be 
simultaneously realized. The optimal choice of the baryon interpolating field seems to be 
around Ioffe's choice. We refer the reader to Refs. [|14] , |T5f for more discussion about the choice 
of baryon interpolating fields. We shall consider the interval —1.15 < t < —0.85 here. For 
t > —0.85 the continuum contributions become large while for t < —1.15 the contributions 
from higher-order terms in the OPE become important relative to the leading-order terms. 

The subscript Sy in Eq. ( |2.ip indicates the presence of the external field. Thus, the 
correlator should be calculated with an additional term 

AC = -S v [u(x)u{x) - d{x)d(x)) , (2.5) 

added to the usual QCD Lagrangian, and — AC added to Hqcb- Since Sy is a scalar 
constant, Lorentz covariance and parity allow one to decompose n(SV, q) into two distinct 



structures 11 



n(sv, q) = n\s v , q 2 ) + u\s v , q 2 )i . (2.6) 

To obtain QCD sum rules, one needs to construct a phenomenological representation for 
n(SV, q) and evaluate n(SV, q) using the OPE. 



A. The dispersion relation and phenomenological spectral ansatz 

To determine the correlator at the hadron level we use the dispersion relation 

ir(W)=f^#<fe (2.7) 

Jo s — q z 

for each invariant function {i = 1, q}, where p*(Sy, s) = -ImrP(Sv> s) is the spectral density. 
Here we have omitted polynomial subtractions which will be eliminated by a subsequent 
Borel transformation. We have also omitted infinitesimal as we are only concerned with 
large and space like q 2 in QCD sum rules. 

In practical applications of QCD sum-rule approach, one usually parametrizes the spec- 
tral density by a simple pole representing the lowest energy baryon state of interest plus a 
continuum which is approximated by a perturbative evaluation of the correlator starting at 
an effective threshold [l|J^,|l^]. When Sy is present, we add — AC to Hqcb, which is equiva- 
lent to increase m u and rrid by Sy and — Sy, respectively. Consequently at the hadron level, 
the baryon spectrum will be shifted. Since we are concerned here with the linear response 
to the external source, Sy can be taken to be arbitrarily small (see below). Thus, there is 
no rearrangement of the spectrum, and we can use a pole plus continuum ansatz for the 
baryon spectral density 
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p\S v , s) = X* B <p l 5(s - M B ) + f?{S y , s)6(s - <) , (2.8) 

where (/)*' = {M B , 1} for {i — 1, g}, and p l (Sy, s) is to be evaluated in perturbation theory. 
Here A B is defined by (OIt/bIB)^ = X B v B with v B the Dirac spinor normalized to v B v B = 
2M B , M B is the mass of the lowest baryon state and Sq is the continuum threshold in the 
presence of the external field. 

Let us now expand both sides of Eq. (|2.7f ) for small Sy 



H* (g 2 ) + S w U\(q 2 ) + ■■■= r ^±ds + S v r ^\ds + ■■■ . (2.9) 

Jo s — q z Jo s — q 

Since Sy is arbitrary, one immediately concludes that 



n (g 2 )= r^- 2 ds, (2.io) 

Jo s — q z 



ni(g 2 )= r^- 2 ds. (2.11) 
Jo s — q 2 

Obviously, Eq. (|2.10| ) leads to the baryon mass sum rules in vacuum which have been ex- 
tensively studied [|1^,|]^,0,[17]] . Here we are interested in Eq. ( j2.11| ), which corresponds to 
the linear response of the correlator to the external source and contains the baryon matrix 
element under consideration (see below). 

Expanding the right-hand side of Eq. (|2.8|) , we find 

Po (s) = A 2 O 5( S - Ml) + pt(s)9(s - 4) , (2.12) 

p\(s) = -2H B M B X 2 B (j)i5'(s - Ml) + AA 2 ^5(s - Ml) 

+A</>* \ 2 B 5(s - Ml) - As* pi{s)5(s - s ) + pi{s)e(s - s ) , (2.13) 

where we have defined 

(2.14) 
(2.15) 
(2.16) 
(2.17) 
(2.18) 

where the first terms are the vacuum spectral parameters in the absence of the external 
field. Note that A0 1 = H B and A0 9 = 0. Treating 5 V as a small parameter, one can use 
the Hellman-Feynman theorem P"B"|,|T9"f to show that 



M* 


= M B + S V H B + • • • 


\* 2 

A B 


= A B + S V AX 2 B + ■ ■ 


si 


= 4 + SyAsl + ■■■ 




= 4> + SyAcf) 1 + ■■■ 


P*\s) 


= f? Q (s) + S y p[(s) + 



(B\uu-dd\B) 
~ 2M B ' (2 ' 19) 
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where we have used covariant normalization (k',B\k,B) = (2Tr) 2 k 5^ 3 \k' — k). 

One notices that p\(s) has specific new features which distinguish it from Pq{s). The 
first term in Eq. ( |2.13| ), which is absent in Pq(s), give rise to a double pole at the baryon 
mass whose residue contains the matrix element of interest. The second and third terms are 
single pole terms; the residue at the single pole contains the information about the transition 
between the ground state baryon and the excited states. In terms of quantum mechanical 
perturbation, the double pole term corresponds to the energy shift while the single pole 
terms result from the response of baryon wave function to the external field. The fouth 
term is due to the response of the continuum threshold to the external source and the last 
term is the continuum contribution. As emphasized in the previous works, the single pole 
contributions are not exponentially damped after the Borel transformation relative to the 
double term and should be retained in a consistent analysis of the sum rules. 

The fourth term has been neglected in Ref. ||11|| . The contribution of this term is sup- 
pressed in comparison with the single pole terms by a factor e _ ^o _M I)/ M2 [see Eqs. ( p.2S - 
|2.36|) 1. If the response of the continuum threshold is small, one can neglect the contribution 
of the fourth term. However, if the response of the continuum threshold is strong, one needs 
to include the fourth term in the calculation. This point has been noticed recently by Ioffe 
in Ref. 



T2|, where a double dispersion relation is considered for the vertex function 



U^q) = / d 4 xe i(I - x (0\T m (x) 



d z(u(z)u(z) — d(z)d(z)) 



%(0)|0) 



(2.20) 



in order to get the appropriate phenomenological representation, 
can be obtained by expanding the right-hand side of Eq. (|2~T 



[This vertex function 
directly] We note that our 



discussion and Eq. ( [2.13|) are consistent with those given in Ref. [0. Substituting Eq. (p. 13 ) 
into Eq. ( [2.11| ), one obtains the appropriate phenomenological representation. 



B. QCD representation 

The QCD representation of the correlator is obtained by applying the OPE to the time- 
ordered product in the correlator. When the external field is present, the up and down quark 
fields satisfy the modified equations of motion: 

(ip -m u - S v )u(x) = , (2.21) 
(ip -m d + Sv)d(x) = , (2.22) 

where p = , y fl (d fl —ig s An) is the covariant derivative. (The equation of motion for the strange 
quark field does not change.) In the framework of the OPE, the external field contributes to 
the correlator in two ways: It couples directly to the quark fields in the baryon interpolating 
fields and it also polarizes the QCD vacuum. Since the external field in the present problem 
is a Lorentz scalar, non-scalar correlators cannot be induced in the QCD vacuum. However 
the external field does modify the condensates already present in the QCD vacuum. To first 
order in Sy, the chiral quark condensates can be written as follows 
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{uu) Sv = (uu) - xSv(uu) , (2.23) 

(dd) Sv = (dd) + X S v (dd) , (2.24) 

(ss)sv = (ss) - XsS v (ss) , (2.25) 

where (O) = (0|O|0). The mixed quark-gluon condensates change in a similar way 

(g s ua ■ Gu) Sv = (g s ucr ■ Qu) - XmS v (g s u(T ■ Qu) , (2.26) 

(g s da ■ Gd)s v = (g s da ■ Qd) + XmS v {g s da ■ gd) , (2.27) 

(g s sa ■ Gs) Sv = (g s scr ■ Qs) Q - x™ s S v (g s so ■ Qs) , (2.28) 

where a • Q = cr^ u Q^ u with Q^ v the gluon field tensor. One can express Xi Xs, Xm, and x ms 
in terms of correlation functions (see Ref. [0). Here we have assumed that the response of 
the up and down quarks is the same, apart from the sign. The Wilson coefficients can be 
calculated following the methods outlined in Ref. |TT|. The results of our calculations for 
the invariant functions U\ and IT^ are given in Appendix A. 



C. Sum rules 



The QCD sum rules are obtained by equating the QCD representation and the phe- 
nomenological representation and applying the Borel transformation. The resulting sum 
rules in the proton case can be expressed as 



- Ci + 3c2 ~ c V m 2 



2H P \ 2 p M 2 - AAp M P M 2 - H p X 2 p M 2 



-M?/M 2 



+ 



AsJmV^ 2 



^2.29) 



i^^aM 4 EnL- 4 / 9 - C4 + C5 - ^ m 2 QM 2 L -26/27 + ^l xa 2 M 2 L 4/9 



4 

ci + 2c 2 
12 



Xm 2 a 



12 

2 r -2/27 Ci — 2c 2 2 2 r -2/27 

L I Xm^Q a L 



12 



2H p \\M P - AAjM 2 



-Ml/Af^ + £i( s g)2 As g M 2 L -8/9 e -4/M^ 

16 



^2.30) 



where a = -4vr 2 (gg) , A 2 = 32tt 4 A£, AA 2 = 32tt 4 AA 2 , and 



(5f s ga • £g) / (gg) . Here 



we have ignored the isospin breaking in the vacuum condensates (i.e., (uOu} ~ (dOd} = 
(qOq) ); the inclusion of the isospin breaking in vacuum condensates only gives small refine- 
ments of the results. We have also defined 
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E 2 = l-e- s °l M * 



+ 



2M 4 M 2 



+ 1 



[2.31] 



and 



C-2 



c 3 = 5r + 2t + 5 , 



c 4 = t 2 + lot + 1 



c 5 



r + 4t + 7 . 



(2.32) 



The anomalous dimensions of the various operators have been taken into account through 
the factor L = ln(M 2 /AQ CD )/ ln(/i 2 /AQ CD ) . We take the renormalization scale \i and 



the QCD scale parameter A QCD to be 500 MeV and 150 MeV [0] 
The sum rules in the S + case are given by 
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3c 2 M 8 E 2 L' 8/9 - 3c 2 xaM % E x + — Xs faM (i E l + ( Cl - 2c 3 )m s aM A E L~ 8/9 

2 



3c 9 



-3c 2 m s faM 4 E L- 8 / 9 + ^m^aM 4 ^- 14 / 27 - ^m s f s m 2 aM 2 L- 38 / 27 

2Ci+3c 2 -6c 3 2_ „,r2 r-38/27 , C l" 2c 3 f/ ,2u2 , 2c 3 2^2 



12 



3 -./Vr.U 2 • ^\m.rr.U^ 



+c 2 xm s /a 2 M 2 + c 2 x s m Ja 2 M 2 

= 2//y A|+M|+ - AAv M s+ M 2 - //v A|+M 2 



-M 2 ,/M 2 



Cl 



-m s (s 



i) 2 L" 4/3 



+ ^ /a^L^ 9 - 3c 2 a4L- 4 / 9 + ^m 2 ^ 26 / 27 A^MV^ 2 > ( 2 . 3 3) 



Scarn.M 6 ^/,- 4 / 3 - ^—^aM^oi -4 '' 9 + 3c 2 faM A E L' A ' 9 - 3c 2X m s aM A E L^ 9 



C3 



4 Xs m s /aM 4 L- 4 / 9 - ^m 2 aM 2 L- 26 / 27 - ^/^aM 2 !" 26 / 27 



5c 2 



+ ^X m m s m 2 aM 2 L~ 26 / 27 - §X ms m J^aM 2 ^ 2 / 27 + ^ X a 2 M 2 L 4 / 9 
-2c 2 x/a 2 M 2 L 4 / 9 - 2c 2 Xs/a 2 M 2 L 4 / 9 - c 2 m s a 2 L- A ' 9 - ^mja 2 !" 4 / 9 



2ci 



C l 2 2 r -2/27 , ^ C 2 /• 2 2 r -2/27 , ^ C 2 /■ 2 2 r -2/27 

~Tn^ m o a L + ^-X/ S m a L 7 + —Xsfm a L ' 



12" 



12 



-^ Xm m 2 a 2 L- 2 / 27 + ^ Xms / s m 2 a 2 L- 2 / 27 + ^ Xm /m 2 a 2 L- 2 / 27 



12 /VI " u 12 

2/A : . A|+M s+ - AAv M 2 



12 



+ 



^(4) 2 - 3c 2 m s a - C jm s fa\As«M 2 L- 8 / 9 e- s »/ M2 , 



(2.34) 



where / = {ss) / {qq) and f s = (g s sa ■ Qs) / (g s q& • Gq)o- The sum rules in the S° case are 
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—M 8 E 2 L- 8/9 + — X aM % E x - 3c 2Xs faM 6 E 1 - 3c 2 m s aM 4 E L- 8/9 
+ ( Cl - 2c 3 )m s faM A E L- 8 / 9 + -^ Xm J W aM 4 E Q L~^ 
-^m s m 2 Q aM 2 L- 38 / 27 - 2c i + g| ~ ^ m j^aM 2 ^ 38 / 27 -^f 2 a 2 M 2 



-c 2 fa 2 M 2 - (d - 2c 3 ) X m s fa 2 M' 2 - ( Cl - 2c 3 ) Xs m s fa 2 M 2 



-3c 2 s faL~ A / g + ^lf s mlaL- 2& l 21 As l e~ s °/ M2 , 



-Ml /M 2 



(2.35) 



3c 2 m s M 6 E 1 L- i/s + jaM 4 E L~ i/9 + 3c 2 faM 4 E L- i/9 - 3c 2X m s aM A E L~ 4/9 
+ 2 -^ Y ^ Xs m s faM*E L-*/ 9 + %nlaM 2 L~ 2 ^ - ? -f f s m 2 aM 2 L~ 2 ^ 
+ ^f Xm m s m 2 aM 2 L- u / 27 + ^ Xms m s f s m 2 M 2 L~ 2 ^ 27 - 2c 2X fa 2 M 2 L^ 9 



-2c 2Xs fa 2 M 2 L" 9 + 2 -^ Xs f 2 a 2 M 2 L" 9 - c 2 m Sj _ 

3 6 



^xfsmVL- 2 ' 27 - ^ Xs ff s m 2 a 2 L- 2 ^ + ^ Xs fm 2 a 2 L- 2 ^ 
+ 5 -§ Xm fm 2 a 2 L- 2 / 27 - ^ Xms ff s m 2 a 2 L~ 2 / 27 + 7 -^ Xms f s m 2 a 2 L^ 27 



f 2 a 2 L-^ 9 - ^-^l ms f s a 2 L^ 9 
5c 2 

7c 2 



2H s o AloM H o - AAlo 



+ 



3c 2 m s a 



-Ml /M 2 



c 3 - 2ci 



m s fa 



As q L- 8 ' 9 M 2 e- s y M2 



(2.36) 



III. SUM-RULE ANALYSIS 



We now analyze the sum rules derived in the previous section and extract the baryon 
matrix elements of interest. Here we follow Ref. [|11] and use only the sum rules Eqs. ( p.30|) , 
( |2.34| ), and ( |2.36| ), which are more stable than the other three sum rules. The pattern 
that one of the sum rules (in each case) works well while the other does not has been seen 
in various external field problems [|^|5T,|T|, |iT||TTf . This may be attributed to the different 
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asymptotic behavior of various sum rules. As emphasized earlier, the phenomenological side 
of the external field sum rules contains single pole terms arising from the transition between 
the ground state and the excited states, whose contribution is not suppressed relative to the 
double pole term and thus contaminates the double pole contribution. The degree of this 
contamination may vary from one sum rule to another. The sum rule with smaller single 
pole contribution works better. We refer the reader to Refs. |7 yT0|JTT| for more discussion 



about the different behavior of various external field sum rules. In the analysis to follow, 
we disregard the sum rule Eqs. ( |2.29| ), ( |2.33Q , and ( |2.35| ), and consider only the results from 
the sum rules Eqs. (|2~3q ), ( ggg ), and ( |2~3E1 ). 

We adopt the numerical optimization procedures used in Refs. fT7| , |2"0|| . The sum rules 
are sampled in the fiducial region of Borel M 2 , where the contributions from the high- 
dimensional condensates remain small and the continuum contribution is controllable. We 
choose 



0.8 < M 2 < 1.4 GeV 2 
1.2 < M 2 < 1.8 GeV 2 



for proton case , 
for S + and H°case 



(3.1) 
(3.2) 



which have been identified as the fiducial region for the baryon mass sum rules |p],|21 
Here we adopt these boundaries as the maximal limits of applicability of the external field 
sum rules. The sum-rule predictions are obtained by minimizing the logarithmic measure 
5(M 2 ) = ln[maximum{LHS,RHS}/minimum{LHS,RHS}] averaged over 150 points evenly 
spaced within the fiducial region of M 2 , where LHS and RHS denote the left- and right-hand 
sides of the sum rules, respectively. 

Note that the vacuum spectral parameters A|, Mb and s l , also appear in the external 
field sum rules Eqs. ( p.29j - f2~30|) and ( |2.33| - f2~ll6| ). Here we use the experimental values for 
the baryon masses and extract A B and s l from baryon mass sum rules using the same 
optimization procedure as described above. We then extract H B , A A 
external field sum rules. 



B , and Asq from the 



0.55 GeV (m u + m d ~ 11.8MeV) [g.[T3| 



m, 



o 



For vacuum condensates, we use a 
0.8 GeV 2 and / ~ /. = 0.8 0,0 

150 MeV The parameter \ has been estimated in Ref. 

perturbation theory gives x — 2.2 GeV -1 . It is also shown that to the lowest order in 5m 
is determined by 



We take the strange quark mass m 



s to be 

The estimate in chiral 

X 



xSm = —7 + 0[(5m) 



(3.3) 



where 7 = (dd) / (uu) — l, and 5m has been determined by Gasser and Leutwyler, 5m/ (m u + 
m^) = 0.28 ± 0.03 p2 |. The value of 7 has been estimated previously in various approaches 



2| with results ranging from — 1 x 10 2 to —2 x 10 3 , which upon using Eq. ( |3.3|) and 



a median value for 5m = 3.3 MeV, corresponds to 

0.5 GeV" 1 < X < 3-0 GeV" 



(3.4) 



We shall consider this range of x values. We follow Ref. |TT| and assume Xm — Xj which is 



equivalent to the assumption that ml is isospin independent. 
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The parameter Xs measures the response of the strange quark condensate to the external 
field, which has not been estimated previously. Since ss is an isospin scalar operator, Xs 
arises from the isospin mixing and we expect Xs < X- Following Ref. flll"| |, one may express 



X s in terms of a correlation function and estimate it in chiral perturbation theory It is 
easy to show that Xs(ss) = ^r(ss} . So, one may determine Xs by evaluating ^r(ss} in 
effective QCD models. Here we shall treat Xs as a free parameter and consider the values of 
X s in the range of < Xs < 3.0 GeV -1 . We also assume that Xms — Xs- 

We first analyze the sum rules for Ioffe's interpolating field (i.e., t — — 1). We start from 
the proton case. The optimized result for H p as function of x is plotted in Fig. [T]. One can 
see that H p varies rapidly with \- Therefore, the sum-rule prediction for the proton matrix 
element H p depends strongly on the response of the up and down quark condensates to the 
external source. (The sum rules in the proton case are independent of Xs and Xsm-) For 
moderate values of x (1-5 GeV -1 < x < 2.0 GeV -1 ), the predictions are 

H p ~ 0.54 - 0.78 . (3.5) 

On the other hand, for large values of x (2.4 GeV -1 < x < 3.0 GeV -1 ), we find H p ~ 
0.97 — 1.25. For small values of x (x < 1.4 GeV" 1 ), the continuum contribution is larger 
than 50%, implying that the continuum contribution is dominant in the Borel region of 
interest and the prediction is not reliable. The predictions for AA^ and Asl also change 
with x in the same way as H p . 

To see how well the sum rule works, we plot the LHS, RHS, and the individual terms 
of RHS of Eq. (|2.30|) as functions of M 2 with x — 1-8 GeV -1 in Fig. |2] using the optimized 
values for H p , A\ p , and Asq. We see that the solid (LHS) and long-dashed (RHS) curves 
are right on top of each other, showing a very good overlap. We also note from Fig. |2| 
that the first term of RHS (curve 1) is larger than the second (curve 2) and third (curve 
3) terms. This shows that the double pole contribution is stronger than the single pole 
contribution and the predictions are thus stable. (Although the second and third terms are 
sizable individually, their sum is small.) 

In Fig. |3|, we have displayed the predicted as function of x for three different values 
of Xs- One notices that is largely insensitive to Xai Du t strongly dependent on x value. 
For x values in the range of 2.2 GeV -1 < x < 3.0 GeV -1 , we find 

Hx+ ~ 1.65 - 2.48 . (3.6) 

For smaller x, we obtain smaller values for Hy,+ . The predictions for AA|+ and Asq change 
in a similar pattern. The sum rule works very well and the continuum contribution is small 
for all x and Xs values considered here. 

The optimized if=o as function of Xs is shown in Fig. [When t — — 1, the sum rule 
Eq. ( |2.36|) is independent of x and Xs-] We see that the result is very sensitive to the Xs 
value. Thus the prediction for if=o has a strong dependence on the response of the strange 
quark condensate to the external field. For moderate Xs (1.7 GeV -1 < 2.2 GeV -1 ), we get 

H s o ~ 1.57- 1.84 . (3.7) 

For larger (smaller) values of Xs, we find larger (smaller) values for H s o. At Xs = 0, we get 
H s o ~ 0.68. The results for AA| and Asl increase (decrease) as Xs increases (decreases). 
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All of the results above use Ioffe's interpolating field (i.e., t = — 1); we now present the 
results for general interpolating field. In Fig. §, we have plotted the predicted H p , H%+, and 
H E o as functions of t for x — 2.5 GeV -1 and Xs = 1-5 GeV" 1 . As t increases, H p , H%+, and 
H s o all increase; the rate of increase is essentially the same for H p and H^,+ , but somewhat 
smaller for if=o. We note that the vacuum spectral parameters Afj and Sq decrease as t 
increases; this leads to a large variation of H p , H%+, and if=o with t. 

The sensitivity of our results to the assumption of Xm = X is displayed in Fig. |6|, where 
t and Xs(— Xms) are fixed at —1 and 1.5 GeV -1 , respectively. The three curves are obtained 
by using x m = X, \Xi and fx, respectively. We note that H p and H s + get larger (smaller) 
as Xm becomes smaller (larger). The results are more sensitive to Xm in the proton case 
than in the S + case. The prediction for H p changes by about 25% while the prediction for 
H%+ changes by about 15% when the Xm value is changed by 50%. This implies that the 
terms proportional to Xm in the sum rules give rise to sizable contributions. The sensitivity 
of our predictions to the assumption of Xsm = Xs is illustrated in Fig. 0, with t — — 1 and 
X — Xm = 2.5 GeV" 1 . The three curves correspond to Xms = Xs, f Xs, and |x s , respectively. 
One can see that both H^,+ and if=o are insensitive to changes in Xms- This indicates that 
the terms proportional to Xms give only small contributions to the sum rules. One also 
notices that depends only weakly on x s - Finally, the effect of ignoring the response of 
continuum threshold is shown in Fig. ^. The solid (dashed) curve is obtained by including 
(omitting) the third term on the RHS of Eq. (|2.30|) . The difference between the two curves 
is large for moderate and large values of x- This shows that the response of the continuum 
threshold can be sizable and should be included in the sum rules. Unfortunately, the response 
of the continuum thresholds has been omitted in all previous works on external field sum 
rules. This was first noticed by Ioffe [12]. 



IV. ESTIMATE OF BARYON ISOSPIN MASS SPLITTINGS 

In this section we estimate the baryon isospin mass splittings using 5m and the baryon 
matrix elements of isovector-scalar current calculated in the previous section. 

The observed hadron isospin mass splittings arise from electromagnetic interaction and 
from the difference between up and down quark masses: 

5m h = (5m h ) e \ + (Sm h ) q , (4.1) 

where (5mh) e \ and (5mfe) q denote the contributions due to electromagnetic interaction and 
due to the up and down quark mass difference, respectively.^] Following Ref. |TT[, one can 



treat 5m as a small parameter and using the Hellman-Feynman theorem to show 

that the octet baryon isospin mass splittings to first order in 5m can be expressed as 



1 This separation is renormalization scale dependent. However, this scale dependence is weak; it 
is thus meaningful to separate the contribution of quark mass difference from that due to electro- 
magnetic interaction (see Ref. |TT||). 
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M n -M p = (M n - M p ) e i + 8mH p , (4.2) 
M s - - M s + = (M E - - M s+ ) el + <5m# s+ , (4.3) 
Af s - - M H o = (M 3 - - M H o) ol + 5mH s o . (A A) 

Note that H n = —H p , H^- = — i?E+j and = —H^o to the lowest order in 8m. There- 
fore, QCD sum rule predictions for H p , H s +, and if=o, along with the electromagnetic 
contributions 



(M n - M p ) el = -0.76 ± 0.30 MeV , (4.5) 
(M s - -M s +) e i = 0.17 ±0.3 MeV , (4.6) 
(M s - - Mso) e i = 0.86 ± 0.30 MeV , (4.7) 

will lead to an estimate of the baryon isospin mass splittings. Taking the experimental mass 
difference [p3]| , one finds 

(M n - M p )^ xp = 2.05 ± 0.30 MeV , (4.8) 
(M s - - M s+ )q xp = 7.9 ± 0.33 MeV , (4.9) 
(Ms- - M H o)^ xp = 5.54 ± 0.67 MeV . (4.10) 

We have seen from last section that the uncertainties in our knowledge of the response of 
the quark condensates to the external field, \ and Xs, leads to uncertainties in the sum- rule 
determination of the baryon matrix elements H-q. (There are also uncertainties in 8m.) 
Therefore, our estimate here are only qualitative. For most of the values for t, x an d Xs 
considered here, the sum-rule analysis gives < H p < H E o < ( see Figs. [5], ^| and 0), 
which implies 

< (M n - M p ) q < (Ms- - M H o) q < (Ms- - M s +) q . (4.11) 

This qualitative feature is compatible with the experimental data. For the baryon interpo- 
lating fields with t — — 1 and moderate x an d Xs values (1.6 GeV -1 < x ^ 2.2 GeV^ 1 and 
1.3 GeV- 1 < Xs < 1.8 GeV- 1 ), we get 

1.95 MeV < (Mn - M p ) q < 2.41 MeV , (4.12) 
4.0 MeV < (Ms- - M s +) q < 6.3 MeV , (4.13) 
4.5 MeV < (Ms- - M H o) q < 5.38 MeV , (4.14) 

where we have used a median value 8m ~ 3.3 MeV. These results are comparable to the 
experimental data, though the result in the £ case is somewhat too small. Smaller and 
larger values of x and Xs lead to correspondingly smaller and larger values for the baryon 
isospin mass differences. As t increases (decreases), the results increase (decrease). 
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V. DISCUSSION 



Our primary goal in the present paper has been to extract the baryon matrix element 
H-q = (B\uu — dd\B) /2Mb for octet baryons. We observe that the sum-rule predictions for 
Hb are quite sensitive to the response of quark condensates to the external isovector-scalar 
field, which is not well determined. This means that our conclusion about H-q can only 
be qualitative at this point. The most concrete conclusion we can draw from this work is 
that QCD sum rules predict positive values for H p , H s +, and if s o and H p < H^o < H^+. 
This qualitative feature is, for the most part, stable against variations of the response of the 
condensates to the external source and the choice of baryon interpolating fields. 

We note that the inequality H p < H^o < indicates SU(3) symmetry violation in 
the baryon matrix elements of the isovector-scalar current. This arises mainly from the 
difference in the baryon interpolating fields used in the QCD sum rules and from the fact 
that the isovector-scalar current is not a SU(3) singlet. Clearly, it is a very interesting topic 
to check this inequality in other effective QCD models. At this stage, it is unclear whether 
the difference in the baryon interpolating fields is connected to the SU(3) symmetry breaking 
in the baryon wave functions. 

In the present study, we derived and used a complete form for the phenomenological 
representation, which has also been given in Ref. This form includes the response 

of the continuum thresholds, which was ignored in Ref. [11]. We found that the neglect 
of the response of the continuum thresholds can have large effect on the extraction of the 
baryon matrix elements. This suggests that the contribution arising from the response of 
the continuum thresholds, neglected in previous works, should be accounted in the study of 
general external field sum rules (see Ref. [^] for estimates of the effects of this contribution 
on the extraction of various physical quantities). 

The spectral parameters in the absence of the external source, M B , A|, and s l , appear 
in all external field sum rules. Unlike the mass, there are no experimental values for the 
coupling A| and the thresholds s l . One usually evaluates these parameters from the mass 
sum rules by fixing the mass at the experimental value. This means that the uncertainties 
associated with the vacuum spectral parameters will give rise to additional uncertainties in 
the determination of the baryon matrix elements of various current, besides the uncertainties 
in the external field sum rules themselves. This is a general drawback of the external field 
sum-rule approach. It is also worth pointing out that it is the product of A| and the baryon 
matrix element appears in the external field sum rules [see Eqs. ( 2.29| — 2.36 )1 • So, it is more 
suitable to determine the product of X B and the baryon matrix element from the external 
field sum rules; one then needs a good knowledge of X B in order to extract the baryon matrix 
element cleanly. 

The sum-rule predictions are fairly sensitive to the choice of baryon interpolating fields. 
This sensitivity arises from both the dependence of the truncated OPE result and the depen- 
dence of the extracted parameters A| and s l on the choice of the baryon interpolating fields. 
We found that the latter has stronger dependence, and hence leads to larger contribution to 
the change of the predictions with t. 

The non-electromagnetic part of the baryon isospin mass difference is essentially given 
by the matrix element Hb multiplied by the light quark mass difference 5m. Given the 
uncertainties in the determination of Hb mentioned above, our estimate of the isospin mass 
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splittings for the octet baryons must be qualitative. It is found that the QCD sum-rule 
predictions yield (M n — M p ) q < (Ms- — M=o) q < (M s ~ — M s +) q . This qualitative result is 
consistent with the experimental data and insensitive to the details of calculation. If we use a 
median value 5m = 3.3 MeV and moderate values for \ an d Xs, we obtain results comparable 
to the experimental values. However, since the response of various condensates to the 
external source and 5m are not precisely known and the uncertainties from other sources 
cannot be accessed systematically, it is not wise to make a critical comparison with data or 
to attempt to extract \ [ an d hence 7 through Eq. (|3.3|)1 and Xs by fitting the experimental 
data. Clearly, further study of the response of the quark condensates to external isovector- 
scalar field is important, along with more accurate determination of the vacuum spectral 
parameters. Effective QCD models may give some independent information on the response 
of the quark condensates while the lattice QCD may offer clean determination of the vacuum 



spectral parameters |L5 



There have been several earlier papers that study the neutron-proton mass difference 
and the baryon isospin mass splittings for other octet baryons p5| , |3lH , based 



on QCD sum-rule approach. In Ref. |[25|| , the baryon mass differences were extracted directly 
from the baryon mass sum rules by including the quark mass difference and the isospin 
breaking in the quark condensates. The contributions of quark-gluon mixed condensates 
were ignored, and somewhat different values for the vacuum condensates and the strange 
quark mass were used. This can lead to large effects on the extraction of the isospin mass 
splittings. The procedure for analyzing the sum rules was also quite different from the 
one used in the present paper. In Ref. @,^], the neutron-proton mass was extracted 
from the difference between the neutron and proton mass sum rules, but the continuum 
contributions were disregarded. In a later calculation |j35| , the authors of Ref. |54j have 



included the continuum contribution in the study of the density dependence of the neutron- 
proton mass difference in the medium. In these works, the contributions from the quark- 
gluon condensates and the change in the continuum thresholds were omitted. The study of 
neutron-proton mass difference in Ref. was based on the mass sum rules directly. Apart 
from keeping the quark mass difference and the quark condensates difference, an attempt 
was made to incorporate the electromagnetic contribution also phenomenologically in the 
sum rules. 



The analysis in Ref. [[Jl[] is more closely related to the present work. The goal of Ref. pl| 
was, however, to determine the parameters 5m and 7 by fitting all isospin mass splittings in 
the baryon octet. The sum rules were obtained for Ioffe's interpolating field by treating the 
quark mass and the isospin breaking in quark condensates as perturbations. On the phe- 
nomenological sides of the sum rules, all spectral parameters, mass, residue, and continuum 
thresholds, were allowed to change. Note that the sum rules derived by us in Sec. || can 
also be derived directly from the mass sum rules. Writing m u = m — 5m/2, m& = m + 5m/ 2 
and assuming \ = —j/Sm [see Eq. ( |3.3| )1, one can differentiate the mass sum rules with 
respect to 5m. For t = —1, one can then identify our sum rules Eqs. ( |2.29H2~30| ) and ( [2.33 - 
|2.36|) with the sum rules given in Ref. |TT]. This coincidence between the sum rules is not 
surprising, since the quark mass term in the QCD Lagrangian can also be regarded as a 
constant external scalar field. We observe, however, that the contributions from dimension 



eight condensates have not been included in Ref. [31]. We have seen in our analysis of the 



sum rules (see Figs. ^| and [I]) that these contributions can be numerically significant. In 
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addition, the authors of Ref . j31| directly used the £ and S mass sum rules from Ref. [16] , 



where all the terms proportional to m u or were neglected. Consequently, some terms 
proportional to m s were not taken into account in the £ and S cases and there was a factor 
two omitted in the contribution from four-quark condensates in the nucleon case. 

We note that the authors of Ref. [31J took a very different procedure in analyzing the 
sum rules. They used both sum rules to eliminate AA| while we used only the more stable 
one. The continuum contribution in sum rules Eqs. ( |2.29|) , ( |2.33| ), and Q2.35Q is large. So, 



these sum rules are likely to be dominated by the single pole terms and the predictions 
based on these sum rules may not be reliable. The size of the continuum contribution was 
not checked in Ref. [[J!]]. Certain assumptions such as Asq = were also used in some 
cases. We notice that the absence of continuum contribution in the external field sum rules 
does not necessarily imply As l = 0. In fact, as long as there is continuum contribution in 
the mass sum rules, one must include As l as a unknown quantity to be determined from 
the sum rules. Any assumption about As l may bypass the information extracted for other 



quantities. The authors of Ref. |3T[ claimed that consistency of the two sum rules can be 
achived for 7 = — (2 ± 1) x 10~ 3 , which is different from the values discussed in the present 
paper (see discussions in Sec. pip . This discrepancy arises mainly from the difference in the 
procedures for analyzing the sum rules. 
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APPENDIX: A 



In this appendix, we give the OPE results for the invariant functions 11} and nf. The 
interpolating fields defined in Eqs. ( |2.2| - [2.4| ) are used in the calculation. We work to the 
leading order in the perturbation theory and to the first order in the strange quark mass 
77V The contributions proportional to up and down quark masses and to gluon condensates 
are neglected. Condensates up to dimension eight are considered. 
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Here ci, C2, C3, C4, and C5 have been denned in Eq. (|2.32|) , and we have ignored the isospin 
breaking in the vacuum condensates (i.e., (uOu) ~ (dOd) () = (qOq) ). All polynomials in 
q 2 , which vanish under the Borel transformation, have been omitted in Eqs. (pO|- |A6|) . 
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FIGURES 



FIG. 1. Optimized sum-rule prediction for H p as function of Xi with Ioffe's interpolating field 
(i.e., t = —1). The other input parameters are described in the text. 



FIG. 2. The left-hand side (solid) and right-hand side (long-dashed) of Eq. ( |2.30| ) as functions 
of Borel M 2 for t = —1, with x = 1-8 GeV -1 and the optimized values for H p , AX p , and Asq. 
The curves 1, 2, and 3 correspond to the first, second, and third terms on the right-hand side of 
Eq. ( pC| ). 

FIG. 3. Optimized sum-rule prediction for as function of x, with t = —1. The three 

curves correspond to Xs = (solid), 1.5 GeV -1 (dashed), and 3.0 GeV -1 (dotted). The other input 
parameters are the same as in Fig. |l[ 

FIG. 4. Optimized sum-rule prediction for H^o as function of Xs> with t = — 1. The other input 
parameters are the same as in Fig. |l[ 

FIG. 5. Optimized sum-rule prediction for H p , H-^+, and H^o as functions of t, with 
X = 2.5 GeV -1 and Xs = 1-5 GeV -1 . The other input parameters are the same as in Fig. [p. 

FIG. 6. Optimized sum-rule prediction for H p and Hj^+ as functions of x, with t = —1 and 
Xs = Xms = 1.5 GeV -1 . The three curves correspond to Xm = X (solid), ^x (dashed), and |x 
(dotted). The other input parameters are the same as in Fig. [I]. 

FIG. 7. Optimized sum-rule prediction for Hj^+ and H-^o as functions of Xs, with t = —1 and 
X = Xm = 2.5 GeV -1 . The three curves correspond to Xms = Xs (solid), |xs (dashed), and |x s 
(dotted). The other input parameters are the same as in Fig. |]. 

FIG. 8. Optimized sum-rule prediction for H p as functions of x> with t = —1. The solid curve is 
obtained by including all three terms on the RHS of Eq. ( 2.30| ) , while the dashed curve is obtained 
by neglecting the third term on the RHS of Eq. ( 2.3C| ). 
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